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INTRODOCTION 

The  use  of  parabolic  equation  (PE)  Bethods 
has  become  very  popular  in  recent  years  tor 
modeling  radar  propagation  effects  in  the 
lower  atmosphere,  especially  for  cases  in 
which  the  vertical  refractive  index  profile 
changes  along  the  propa'^ation  path.  The  PE 
method  normally  used  is  the  split-step  method 
described  by  Tappert  (1),  which  has  been 
implemented  by  Dockery  (2),  Craig  and  Levy 
(3),  and  others.  An  advantage  of  the  PE 
method  is  its  ability  to  compute  propagation 
effects  within  the  horizon  as  well  as  beyond 
the  horizon,  thereby  allowing  computations  to 
be  made  in  all  regions  of  practical  interest 
to  radar  engineers  or  operators  with  just  one 
model.  However,  a  significant  disadvantage  of 
the  split-step  PE  method  is  that  it  requires 
extensive  computation.  Computational  reejuire- 
ments  increase  with  higher  frequencies, 
larger  antenna  beamwldths,  and  higher 
altitudes  for  which  results  are  desired.  For 
many  practical  combinations  of  these 
parameters,  the  use  of  PE  models  on  personal 
computers  is  impractical  without  extra 
hardware  such  as  transputers. 

This  paper  presents  a  hybrid  propagation 
model  called  the  Radio  Physical  Optics  <RPO) 
model  that  uses  a  combination  of  ray  optics 
(RO)  and  split-step  PE  methods  to  overcome 
the  high  computational  burden  of  pure  split- 
step  PE  methods.  HPO  considers  four  regions 
shown  in  Figure  1,  At  ranges  less  than  2500 
meters  and  for  all  elevation  angles  above  5 
degrees,  RPO  uses  a  flat  earth  (FE)  model 
that  ignores  refractive  and  earth-curvature 
effects.  For  the  region  beyond  the  FE  region 
where  the  grazing  angles  of  reflected  rays 
from  the  transmitter  are  above  a  small 
limiting  value,  a  full  RO  model  is  used  that 
accounts  for  the  effects  of  refraction  and 
earth  curvature.  The  PE  model  is  used  for 
ranges  beyond  the  RO  region,  but  only  for 
altitudes  below  a  maximum  PE  altitude  deter¬ 
mined  by  the  maximum  1024-point  fast-Fourier 
transform  (FFT)  allowed.  For  ranges  beyond 
the  RO  region  and  heights  above  the  PE  re¬ 
gion,  an  extended-optics  (XO)  method  is  used 
that  is  initialized  by  the  PE  model  at  the 
maxiaua  PE  altitude,  a.nd  uses  ray-optics 
methods  to  propagate  the  signal  to  higher 
altitudes.  Continuity  of  the  solutions  across 
each  region's  boundaries  is  kept  to  less  than 
0.1  dB  by  careful  selection  of  the  limiting 
RO  grazing  angle  and  the  maximum  PE 
propagation  angle. 

Figure  2  shows  an  RPO  example  for  3300  MHz 
with  a  transmitter  at  30.5  meters  above  sea 
level.  The  figure  shows  propagation  loss  In 
dB  on  a  gray-shade  scale  versus  range  and 
height.  Propagation  loss  is  the  ratio  of 
transmitted  to  received  power  using  the 
actual  transmitter  antenna  pattern  but 
normalized  to  0  dB  antenna  gain.  The 
environment  is  a  measured  range-dependent 
case  from  Point  Loma  in  San  Diego,  California 
on  a  path  toward  Cu.adalupe  Isle,  Mexico.  The 


refractivity  profiles  along  this  path 
indicated  a  trapping  layer  that  increase-.!  i 
height  with  increasing  range,  such  tliat  a 
surface-based  duct  existed  at  Point  Ixima  and 
a  low-elevated  duct  existed  at  250  km.  As  a 
reference  point  to  Figure  1,  ttie  love;  left 
corner  of  the  XO  region  in  this  example  is  at 
a  range  of  95  km  a.nd  a  height  of  1260  m. 

a!ix_Am,csusopEi, 

Refractive  index  profiles  are  specified  in 
RPO  at  one  or  more  ranges,  and  are  entered  in 
terms  of  modified  refractivity,  H,  defined  as 

M  •  (n-1 -z/a)  xlO-*  ll! 

where  n  is  the  radio  refractive  index,  z  is 
height,  and  a  is  the  earth's  radius. 
Computations  in  the  RO  region  do  not  consider 
range-dependent  environmental  effects,  but 
rather  depend  only  on  the  vertical  M  profile 
at  the  transmitter. 

The  key  to  RPO's  efficiency  is  keeping  the 
angles  considered  in  the  PE  model  as  small  as 
possible  by  maximizing  the  RO  region.  A 
limiting  grazing  angle  Rj  for  reflected  rays 
is  determined  that  defines  the  maximum  range 
and  altitude  to  which  the  RO  method  can  be 
applied.  R,  is  first  computed  as 

I)!,  •  0.04443/f''’  ti) 


where  is  in  radians  and  f  is  frequency  in 
HHz.  R,  from  equation  (2)  is  2.5  times  the 
limit  given  by  Reed  and  Russell  (4).  The 
factor  2.5  was  chosen  to  ensure  that  ctrort. 
in  the  RO  solution  would  be  less  than  O.I  dP, 
♦„  from  equation  (2)  is  limited  to  values 
above  0.002  radians,  and  then  doubled  if  a 
range-dependent  environment  has  been 
specified.  Finally,  Rj  is  increased  by  an 
amount  SR  to  account  for  ducting  given  by 

S*  .  y/'2xlO'‘(H,-M.7 

where  is  M  at  the  surface  and  in  the* 
minimua  value  of  M  at  all  heights. 

The  RO  nethod  consists  of  tracing  a  senes  of 
direct  and  reflected  rays  through  selected 
control  points,  and  then  interpolating  the 
magnitudes  of  the  direct  and  reflected  rays 
and  the  phase  angle  between  then  at  each 
desired  RPO  output  point.  The  magnitude  of 
each  ray  is  computed  from  a  spreading  term 
relative  to  free-space  spreading,  and  the 
phase  angle  is  determined  from  the  optical 
path  length  differences  from  the  ground  range 
for  each  ray.  The  direct  and  reflected  ray 
through  a  given  point  are  characterized  by 
their  elevation  angles  at  the  transmitter,  f, 
and  o  ,  respectively.  The  raytrace  method  ir; 
summarized  below  for  a  single  rtep. 

The  elevation  angle  at  the  end  of  the  step, 
a^,  is  given  by 


f 


(4) 


Where  o,  is  the  elevation  angle  at  the 
beginning  of  the  step,  and  and  M,., 
represent  H  at  the  beginning  and  end  of  the 
step.  The  range  increment  Ax,  spreading 
increment  AS,  and  optical  path  length 
difference  increment  AD  over  the  step  are 
given  as 

Ax  »  (a,-o,),'g, 


AS  =.  (o/a,-a/ao)/gi  (5) 

ad  -  [  {10-^,-aJ/2)  (a,-ao)  •  (aJ-o’)/3  )/q, 


where  a  is  the  elevation  angle  at  the 
transmitter,  and  z,  and  z,,,  are  the  heights 
at  the  beginning  and  end  of  the  step.  The 
total  range,  x,  total  spreading  term,  S,  and 
total  optical  path  length  difference,  D,  are 
given  by  the  sums  of  Ax,  AS,  and  AD  over  all 
steps  along  each  ray.  The  propagation 
factors,  and  F^,  for  the  direct  and 
reflected  rays  and  the  total  phase  angle 
between  the  rays  in  radians,  n,  are  given  by 

^  ^  I 

^0 


F^ 


ff 


n  -  (D,-D,,)  k  *  <t> 


(*) 


where  x  and  z  are  the  t.ii'../e  and  heig.ht;  rj,;,.';, 
the  transmitter  and  is  the  t  ransaii  t  tei 
height  above  the  surface.  The  g!,izing  angle, 
♦  ,  and  the  direct  ray  elevation  angle  at  trie 
transmitter,  a^,  are  given  by 


♦  ■*  tarr  ' !  (  z  ,  X  ! 
-  tan  '( 


(  10  i 


and  the  reflected  ray  eievation  a.ngle  is 
given  by  a,  -  -  $.  The  total  piurse  lag,  n. 
In  radians  is 


n  ‘  (Ti  r,}  y.  ■  .»  .  (13) 

The  propagation  factor,  Y,  tar  the  coherent 
sum  of  the  two  ray  components  is  computed  by 
equation  (7),  where  F^  ■*  f j  and  Y\  -  f,B. 
Propagation  loss  is  computed  from  equation 
(8),  substituting  r,  for  x. 


The  split-step  PE  model  follows  Doclcery  (2), 
where  the  complex  field  u(x,z)  is  advanced  to 
U(X4-4X,Z)  by 

ru  J 

u(x*4x,  z)  '  (U  (X.  p)  e  ) 

where  M  is  the  modified  refractivlty  as 
defined  in  equation  (1).  The  Fourier 
transform  /  of  u(x,z)  is  defined  as 

U{x,p)  -  .?'(u(x,  z)  1  ■  J  ’  u(x,  z)e 


where  and  f,  are  antenna  pattern  factors 
corresponding  to  the  elevation  angles  at  the 
transmitter,  and  a,,  respectively.  0^  and 
0^  are  elevation  angles  at  the  terminal  point 
of  each  ray  for  the  direct  and  reflected 
rays,  k  is  the  wave  number  2w/X  where  A  is 
wavelength  in  meters.  R  and  <(>  are  the 
magnitude  and  phase  lag  of  the  reflection 
coefficient,  which  are  computed  for 
horizontal,  vertical,  or  circular 
polarization  and  include  the  effect  of  sea 
roughness  based  on  wind  speed  in  the  same 
manner  as  described  by  Patterson  et  al.  (5). 
The  sum  of  the  two  ray  components  is  given  by 

F^  -  F^  ♦  fJ  -  2F„F,  cos  n  (f) 


where  F  is  the  propagation  factor  defined  as 
the  ratio  of  the  field  strength  to  the  free- 
space  field  strength.  Propagation  loss  L  in 
dB  is  computed  as 

<8) 

L  =■  20  log  f  *  20  log  X  -  10  log  F^  -  27.56 

where  log  is  base  10,  f  is  frequency  in  MHz, 
and  X  is  ranqe  in  meters, 

FIAT  BRRTH  MODEL 

The  flat  earth  model  ignores  all  effects  from 
refraction  and  earth  curvature,  and  computes 
the  direct  and  reflected  ray  interference 
pattern  using  straight  line  paths.  Full 
account  is  given  to  the  antenna  pattern  and 
sea-surface  reflection  coefficient,  including 
rough  surface  effects.  The  direct  path 
length,  r^,  and  the  reflected  path  length, 

Tj,  are  given  by 

r,  = 

<») 


where  p  “  k  sin  9,  and  0  is  the  angle  from 
the  horizontal.  In  RPO,  only  real-valued  sine 
FFTs  are  used,  with  which  the  real  and 
imaginary  parts  of  u  are  transformed 
separately.  A  filter  is  applied  to  the  upper 
1/4  of  the  field  in  both  z-  and  p-spaces  at 
each  step  to  ensure  that  the  field  reduces  to 
zero  at  the  top  of  the  transform.  Transform 
size  varies,  but  never  exceeds  1024  points. 

The  starting  field  at  x  *  o  is  constructed  in 
p-space  based  on  image  theory  and  far-field 
approximations.  Thus 

0(0, p)  •  G  [f,e  -  f.Re""'-*’] 

which  is  normalized  by  G  such  that 
F^  "X  ju(x,z)  1' 

which  is  used  in  equation  (8)  to  compute 
propagation  loss  in  dB.  The  magnitude  R  and 
phase  lag  <p  oC  the  reflection  coefficient 
include  polarization  and  rough  surface 
effects  in  the  same  manner  as  used  in  the  BO 
region.  A  Gaussian  taper  function,  much 
stronger  than  the  filter  referred  to  above, 
is  applied  to  the  upper  1/4  of  the  starting 
solution.  T.tis  fu.-.cti^n  reduces  the  field  up 
to  70  dB,  in  addition  to  any  reduction  by  the 
antenna  pattern. 

The  selection  of  the  height  increment  6z 
between  adjacent  values  of  u  and  the  range 
increment  4x  follows  methods  outlined  by 
Tappert  (1),  and  summarized  an 

4  z  -  In, 


4  X  "  2  V.  (  4  ' 


-  J(Z»Z,)'*X^ 


(16) 


-where  is  an  appropriate  integer  nuinijer, 
and  6^^  is  the  aaximuo  angle  that  can  be 
accoauaodated  in  the  PE  solution.  In  RPO, 
is  detenained  as  the  maxlmua  Integer  such 
that  3/4  of  8  just  exceeds  the  oiaximuiB 
elevation  angle  of  the  RO  limiting  ray  at  any 
height  below  3/4  of  the  maximum  PE  height. 

eXTENDEO  OPTICS  HODEt 

The  propagation  factor  F  is  computed  from  the 
PE  solution  at  the  top  of  the  PE  region  and 
is  used  to  initialize  an  RO  model  In  the  XO 
region.  This  model  traces  rays,  along  which  F 
from  the  PE  solution  is  held  constant,  taking 
account  of  the  full  range-dependent 
refractivity  environment,  and  based  on 
initial  angles  ^  at  the  PE/XO  boundary.  If  a 
reflected  ray  from  the  transmitter  exists  at 
the  PE/XO  boundary,  then  0  is  the  local 
elevation  angle  of  this  ray  at  the  PE/XO 
boundary.  The  greatest  range  at  which 
reflected  rays  exist  along  the  PE/XO  boundary 
is  called  the  optical  limit.  For  ranges 
beyond  the  optical  limit,  is  given  by 

^  =■  •  2><10-‘(M  -  M^) 

whera  0^  and  are  0  and  modified 
refractivity  at  the  optical  limit,  and  M  is 
the  modified  refractivity  at  the  desired 
range.  Using  linear  interpolation  techniques, 
F  is  defined  at  all  points  within  the  XO 
region,  and  propagation  loss  is  calculated 
using  equation  (8)  . 

COMPARTSOW  TO  OTHER  METHODS 

Figure  3  compares  RPO  results  to  results  from 
a  waveguide  program  and  a  split-step  PE 
model.  The  case  presented  is  a  homogeneous 
surface  based  duct  defined  by  the  modified 
refractivity  profile  of  Table  1.  The 
frequency  is  3000  MHr ,  the  antenna  height  is 
30.5  meters,  the  polarization  is  horizontal, 
the  antenna  pattern  is  omni-directional,  and 
a  smooth  sea  surface  is  assumed.  Figure  3 
plots  propagation  factor  in  dB  versus 
altitude  at  a  range  of  185  km.  Figure  3a 
shows  the  RPO  results  for  the  ducting  case  by 
a  solid  curve,  and  the  corresponding  standard 
atmosphere  case  by  a  dot-dash  curve.  The 
dotted  horizontal  lines  in  Figure  3a  indicate 
the  boundaries  of  the  PE,  XO,  and  RO  regions 
for  the  ducting  case  at  the  185  km  range.  The 
effects  of  this  strong  duct  in  all  three  of 
the  RPO  regions  is  clear  from  a  comparison  of 
the  ducting  and  standard  curves.  Figure  3b 
presents  waveguide  results  for  this  same  case 
previously  reported  by  Hitney  et  al.  (6). 

The  waveguide  results  are  practically 
ide.ntical  to  the  RPO  results.  Figure  3c 
presents  results  from  a  split-step  PE  model 
described  by  Barrios  (7).  These  results  are 
also  virtually  identical  to  the  RPO  results. 

The  times  required  to  compute  the  ducting 
cases  of  Figure  3  on  a  25  MHz  IDM/PC- 
compatible  486  computer  were  69,  381,  and  310 
seconds  for  the  RPO,  waveguide,  and  PE 
models,  respectively.  In  more  stressful 
cases,  RPO  has  proven  to  be  25  to  100  times 
faster  than  split-step  PE  models.  For  some 
practical  combinations  of  higher  heights  and 
frequencies  and  using  beamwidths  above  a  few 
degrees,  PE  models  are  usually  impractical 
due  to  large  transform  size  requirements.  Rfv' 
has  no  difficulty  with  these  applications. 

One  limitation  of  the  current  implementation 
of  Pl'O  IS  the  use  of  a  range-independent  RO 
model  in  range-dependent  envi ronmentr. ,  which 


in  extreme  cases  can  cause  dis;::;,:'.! 
along  the  RO/PE  boundary.  However,  expes lere u 
has  shown  that  for  realistic  env  i  r  s  , 

such  as  the  case  of  Figure  2,  tire  limit  in  j 
grazing  angle  is  high  enough  to  ensure 
acceptable  RO  results.  A  second  1  i  si  t.at  ic-n 
concerns  rough  surface  effects  in  the  PE 
model,  which  are  included  in  the  st.ittinj 
solution  to  account  for  higher  .ingle  effe.  tn 
but  are  ignored  thereafter.  Comparisons  of 
waveguide  and  RfO  result-s  indic.ite  t  fi  i  s 
simplification  can  affect  some  higii  liegueri,, 
applications  in  ducting  environments 

cetietcsjteiis 

Based  on  comparisons  of  results,  R!wd  Is  as 
accurate  as  waveguide  or  split-step  PE 
models.  RPO  can  compute  results  for  range- 
dependent  environments  In  the  same  way  PE 
models  can,  but  is  applicable  to  wider 
antenna  patterns  than  typical  PE  Bodels.  RKj 
has  proven  to  be  much  faster  than  split-step 
PE  models,  yet  it  requires  far  less  memory 
than  typical  PE  models.  The  accurate  results, 
wider  applicability,  faster  computation 
times,  and  smaller  memory  requirements  make 
RPO  an  Ideal  model  for  use  in  all  propagation 
assessment  or  engineering-aid  programs, 
including  those  hosted  on  personal  computers. 

Special  thanks  are  due  to  F.D.  Tappert  for 
providing  a  compact  working  PE  model  in 
August  1990,  from  which  RPO's  PE  model  was 
derived.  The  development  of  RPO  was  sponsore-,! 
by  the  office  of  Haval  Technology.  Patent 
protection  Is  pending  on  this  work.  This 
paper  Is  in  the  public  domain. 
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Tabl*  1.  Modified  re/ractivity  versus  height 
profile  for  the  ducting  case  of  Figure  3. 


